NASA TECHNICAL NOTE 


LO 



NASA IN D-6075 

C,/ 




LOAN COPY: RETURN !^j 


AFWL (DOGL) 
KIRTLANP AFB, N- ^ 


O" I 
tr ! 


ANALYSIS AND COMPUTER PROGRAM 
FOR EVALUATION OF ROTATING 
INCOMPRESSIBLY LUBRICATED 
PRESSURIZED THRUST BEARINGS 

by David P. Fleming 

Lewis Research Center 
Cleveland^ Ohio 44133 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION • WASHINGTON, D. C. « FEBRUARY 1971 


I 


TECH LIBRARY KAFB, NM 



TECH LIBRARY KAFB, NM 



1. Report No. 2. Government Accession No. 

NASA TN D-6075 

3. Recipient's Catalog No. 

4. Title and Subtitle ANALYSIS AND COMPUTER PROGRAM FOR 
EVALUATION OF ROTATING ESTCOMPRESSIBLY LUBRICATED 
PRESSURIZED THRUST BEARINGS 

5. Report Date 

February 1971 

6. Performing Organization Code 

7. Author(s) 

David P. Fleming 

8. Performing Organization Report No. 

E-5938 

10. Work Unit No. 

129-03 

9. Performing Organization Name and Address 

Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio 44135 

11. Contract or Grant No. 

13. Type of Report and Period Covered 

Technical Note 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D. C. 20546 

14. Sponsoring Agency Code 

1 15. Supplementary Notes I 


16. Abstract 

An analysis and FORTRAN IV computer program are presented to enable rapid evaluation of pres- 
surized thrust bearing designs using an incompressible lubricant. Included in the analysis are the 
effects of two self-acting journal bearings which may be used to provide a radial load capacity. 
Bearing load, torque, lubricant flow rate, and other quantities of interest are calculated. Either 
orifice or capillary restrictors may be used and effects of bearing rotation are included. Program 
input and output can be in U. S. Customary or metric (SI) units. Analytical predictions agreed 
well with experimental data from a series -hybrid fluid -film rolling-element bearing. 


17. Key Words (Suggested by Author (s) ) 18. Distribution Statement 

Bearings Hydrostatic bearings Unclassified - unlimited 

Externally pressurized Thrust bearings 

bearings 


1 

19. Security Classif. (of this report) 

Unclassified 

20. Security Classif. (of this page) 
Unclassified 

21. No. of Pages 
28 

22. Price* 

$3.00 


For sale by the National Technical Information Service, Springfield, Virginia 22151 


ANALYSIS AND COMPUTER PROGRAM FOR EVALUATION OF ROTATING 


INCOMPRESSIBLY LUBRICATED PRESSURIZED THRUST BEARINGS 

by David P. Fleming 
Lewis Research Center 

SUMMARY 

An analysis and FORTRAN IV computer program are presented to enable rapid eval- 
uation of pressurized thrust bearing designs using an incompressible lubricant. Included 
in the analysis are the effects of two self-acting journal bearings which may be used to 
provide a radial load capacity. Bearing load, torque, lubricant flow rate, and other 
quantities of interest are calculated. Either orifice or capillary restrictors may be used 
and effects of bearing rotation are included. Program input and output can be in U. S. 
Customary or metric (SI) units. Analytical predictions agreed well with experimental 
data from a series-hybrid fluid-film rolling-element bearing. 


INTRODUCTION 

For the evaluation of a series -hybrid fluid -film ballbearing (ref. 1), the investigators 
used, for the fluid-film bearing, a pressurized thrust bearing in conjunction with a self- 
acting journal bearing. The lubricant was the same as that used for the ball bearing, 
namely, a type II ester oil. In the series -hybrid bearing (fig. 1), a fluid -film bearing 
and a ball bearing are coupled in series. Each component bearing carries the full system 
load, but the two bearings share the speed. One element of the fluid-film bearing rotates 
at shaft speed. The other element of the fluid-film bearing rotates with the inner race of 
the ball bearing at a speed less than shaft speed. The ball bearing outer race is station- 
ary. The intermediate member rotates at a speed such that the torques of the fluid-film 
and ball bearings are equal. As shown in reference 1, the series -hybrid bearing has the 
potential of substantially increasing rolling-element bearing fatigue life at high speed. 

Since the speed sharing between the fluid-film and ball bearings depends on the torque 
characteristics of the two component bearings, the bearing designer must properly size 
the fluid -film bearing to get a useful reduction in ball bearing speed and provide ade- 
quate load capacity. 

Self-acting journal bearings have been quite thoroughly investigated (see, for exam- 
ple, refs. 2 and 3). Surprisingly, there appears to be no published information on rotat- 
ing, compensated, pressurized thrust bearings using incompressible lubricants. Refer- 



ence 4 gives design information for hydrostatic bearings with relative motion between the 
bearing parts. With liquid lubricants, however, centrifugal effects due to bearing rota- 
tion can significantly alter the bearing’s characteristics. Formulas have been given for 
the load capacity of rotating uncompensated pressurized thrust bearings (for example, 
ref. 2, p. 203). The disadvantage of this type of bearing is that it has no stiffness and 
therefore cannot accept varying loads. It is, fortunately, a straightforward matter to 
add the effects of compensating restrictors (orifices or capillaries) to the analysis of an 
uncompensated pressurized bearing. 

The objectives of this report are the following: 

(1) To analyze an orifice (or capillary) compensated thrust bearing, including the ro- 

tational effects and the effects of any journal bearings adjoining the thrust bearing 

(2) To present a digital computer program to carry out the analysis 

(3) To give sample results of the analysis, and to compare these results with experi- 

mental data from the series -hybrid bearing 

SYMBOLS 

Cj orifice discharge coefficient 

d 

Cj radial clearance in inner journal bearing, in. (m) 

Cq radial clearance in outer journal bearing, in. (m) 

d orifice or capillary diameter, in. (m) 

h thrust bearing clearance, in. (m) 

L journal bearing length, in. (m) 

Lp capillary length, in. (m) 

n number of restrictors 

2 2 

p pressure, Ib/in. (N/m ) 

Q lubricant flow rate, in. /sec (m /sec) 

R radius, in. (m) 

Re film rotational Reynolds number, r|wj - Wglhp/p 
r radial coordinate, in. (m) 

T bearing torque, in. -lb (N-m) 

V velocity, in. /sec (m/sec) 

W bearing load, lb (N) 


2 



z axial coordinate, in. (m) 

e journal bearing eccentricity ratio 

9 circumferential coordinate, rad 

2 2 

II lubricant dynamic viscosity, Ib-sec/in. (N-sec/m ) 

O O 

p lubricant density, Ib/in. (kg/m ) 

u) angular velocity, rpm or rad/sec 

cOq mean angular velocity, eq. (7), rpm or rad/sec 

Subscripts: 

c restrictor exit 

i inner 

j journal bearing 

min minimum 

0 outer 

p pocket 

r radial direction 

s supply 

z axial direction 

& circumferential direction 

1 upper thrust surface 

2 lower thrust surface 


ANALYSIS 

The bearing to be analyzed appears in figure 2. It comprises a thrust bearing and 
two journal bearings at the inside and outside thrust bearing radii. The journal bearings 
enter into the analysis only as they contribute additional constant resistances to the 
throughflow of lubricant and as they add to the bearing torque. Their lengths may be set 
to zero when they are not present. A circle of orifices is at radius R^. Alternatively, 
the bearing may have capillary restrictors at radius R^. The number of orifices or cap- 
illaries is assumed large enough to constitute a line source of lubricant. The thrust face 
clearance varies with radius in a stepwise manner as shown. The lubricant is supplied 
to the orifices at pressure p^; it leaves the bearing at the reference pressure p = 0. 
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The starting point of the analysis is the Navier -Stokes equations for incompressible flow 
(ref. 5). The following assumptions are then imposed: 

(1) The flow is laminar and steady with time. 

(2) Rotational symmetry prevails, that is, 9/30 = 0. 

(3) Pressure in the thrust bearing is constant in the axial direction. 

(4) There are no body forces acting on the fluid. (Centrifugal force appears as an 
inertia term. ) 

(5) The axial velocity is negligibly small in the thrust bearing. 

With these assumptions, the Navier-Stokes equations become, in cylindrical coordinates. 
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The continuity equation has been incorporated into equation (la). In a rotating bearing, 
the radial fluid velocity v^. is much less than the circumferential velocity v^. Thus 
first term of equation (la) can be neglected relative to the second. It is assumed that 
viscous forces predominate over inertia forces. Thus the first term of equation (lb) can 
be eliminated. (These are standard assumptions in lubrication work. See, for example, 
ref. 3, p. 68). Equations (la) and (lb) become 



dr 


M 



9r 




A solution of equation (2b) is 


Mq = rw2 - r(o)2 - coj) | 


(2a) 


(2b) 


( 3 ) 
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where z is measured from the lower thrust surface (fig. 2). Equation (3) may be sub- 
stituted into equation (2a) and the result integrated twice in z. After the boundary con- 
ditions Vj. = 0 at z = 0 and = 0 at z = h are applied, the result is 


J- ^ . pr 

2ju dr p 
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Lubricant Flow Rates 


The quantity of lubricant flowing radially at any radius r is found from 


•/' 

*'0 


Q = 27rr / v^dz 


Substituting equation (4) into equation (5) and integrating yield 


Q = 


Trrh'^ 




where is defined by 
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(5) 


( 6 ) 


(7) 


The speed o.^, which may be regarded as an average of oJj and 0 - 2 , is used to calculate 
rotational effects. Equation (6) may now be integrated with respect to r to find the relation 
between flow and pressure. The clearance h has been assumed to vary in a stepwise man- 
ner. Thus, integration over any interval of constant clearance is straightforward. For 
example, integrating from r = R to r = R R gives 
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Pressure drop through the journal bearings is easily calculated from the expression 
for flow in a narrow slot (ref. 3, p. 99). For the outer journal bearing, neglecting the 
effect of eccentricity, 


^o ~ ^ref 


7TR 

o o 


( 9 ) 


With this, total lubricant flow through the bearing can now be given as that flowing 
inward from the orifices and that flowing outward: 



The negative of Qj is taken because, at any radius, outward flow is defined as positive. 

The total flow Q must now be matched to the flow through the restiictors. For ori- 
fice restrictors, this is (ref. 4, p. 103) 


Q = C , nff — 
^ 4 


2(Ps - Pc) 


(11a) 


and for capillary restrictors , 


n7Td‘^(Pg - Pp) 
128 juLj^ 


(11b) 


Equation (11a) or (11b), as appropriate, may be combined with equation (10), and the 
pressure p^ downstream of the restrictors solved for algebraically. The flow rate Q 
can then be found from equation (10). 
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Bearing Loads 


The thrust bearing load is given by 


W = 



pr dr 


( 12 ) 


The pressure p is given by equation (8) for r R^, and by similar expressions 

for other radii. These expressions may be substituted into equation (12) to yield the 
load. For example, 



Inner and outer journal bearings (fig. 2) may be used to provide a radial load capacity. 
These are assumed to be purely self-acting bearings. For small length-to-diameter 
ratios, the load is adequately given by the short bearing approximation. From reference 
3, (p. 84), for the inner bearing. 


W.. = 
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MRj 


^1 - ^^2 


4Cf 






(- 






+ 16ef 


2 ] 1/2 


(14) 


where is the journal -bearing eccentricity ratio. 


Bearing Torque 


The velocity Vq varies linearly across the film (eq. (3)). Thus, thrust bearing 
torque is easily calculated. For the annulus from r = R to r = R„„, 



(15) 
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Similar expressions apply to the other sections. For the inner journal bearing, neglec- 
ting the effect of eccentricity on torque, the torque is 

^ _2^R?LiPlcUj - CU 2 I 

c ‘ 

'"i 

A similar expression applies to the outer journal bearing. Total torque is merely the 
sum of the various component torques: 


T = + Tp, + Tj + 


T. + T.. 
]o ]i 


(17) 


COMPUTER PROGRAM 

Though the preceding analysis is simple and straightforward, a considerable effort 
would be needed to calculate the large number of numerical results needed in evaluating 
several bearing designs and operating conditions. To reduce this effort, a computer pro- 
gram was written. The program is written in FORTRAN IV for the IBM 7094-11 com- 
puters at Lewis Research Center. It can be used on any digital computer accepting 
FORTRAN IV programs. The program listing appears in the appendix. Figure 3 is a 
flow chart of the program. Four options are available in the program: 

(1) Input and output may be in either U. S. Customary or metric (SI) units. To use 
metric units, set METRIC = .TRUE, in the input list. 

(2) Orifice or capillary restrictors may be used. If the capillary length LR is not 
zero, capillary restrictors are assumed; if LR = O, orifice restrictors are assumed. 

(3) Output may be for a specified set of clearances or a specified set of bearing 
torques. If TMAXX is zero, output will be for NH outer thrust face clearances (h^) 
starting at h^ = HOO and proceeding up by increments of DH. If TMAXX is nonzero, 
the program will still start with h^ = HOO. It will then increase h^ by increments of 
DH until the torque decreases to TMAXX. (HOO must be chosen so the corresponding 
torque is greater than TMAXX). It will then make calculations for clearances such that 
bearing torque drops by increments of DT from TMAXX to TMIN. Under either op- 
tion, the program stops incrementing clearance or torque whenever the bearing load be- 
comes negative. 

(4) Supply pressure PS may be specified. Or, if PS is not specified, the pro- 
gram will calculate a supply pressure based on the pressure available from centrifugal 
force at the radius of the restrictors: 
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( 18 ) 




To have the program calculate p , set CENT = .TRUE. 

s 

As a check on whether the assumed laminar conditions prevail, the program calcu- 
lates Reynolds numbers based on the thrust bearing film thickness and rotational speed: 

Re = r] Wj - (19) 

The three Reynolds numbers calculated are for the outer thrust face land (using R^ and 
h^), for the inner thrust face land (using R^. and hj), and for the pocket (using R^^ 

and hp). If any one of these Reynolds numbers exceeds 1000, the flow is probably tur- 
bulent in that region (ref. 6). Computed results will then be inaccurate. 


Program Input 

All input is in NAMELIST format using the NAMELIST name BRG. Sample input ap- 
pears in figure 4. Note that all data cards begin in column 2. Following is a description 
of the input variables: 


FORTRAN 

Analysis 

Description 

name 

symbol 


CD 

Cd 

orifice discharge coefficient 

CENT 


if . TRUE. , PS calculated from eq. (18) 

Cl 

Ci 

radial clearance in inner journal bearing, in. (m) 

CO 

Co 

radial clearance in outer journal bearing, in. (m) 

D 

d 

diameter of orifices or capillaries, in. (m) 

DH 


increment in thrust bearing clearance, in. (m) 

DT 


increment in bearing torque, in. -lb (N-m) 

DW2 


increment in speed 0 ^ 2 , rpm 

HII 


initial value of inner thrust clearance h., in. (m) 

HOO 


initial value of outer thrust clearance h^, in. (m) 

HPP 


initial value of pocket thrust clearance hp, in. (m) 

LI 

Li 

length of inner journal bearing, in. (m) (may be zero) 
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FORTRAN 

Analysis 

Description 

name 

symbol 


LO 

Lo 

length of outer journal bearing, in. (m) (may be zero) 

LR 

^R 

length of capillary restrictors, in. (m) (zero for orifice 
restrictors) 

METRIC 


if . TRUE. , input and output are in metric (SI) units 

MU 


0 0 

lubricant dynamic viscosity, Ib-sec/in. (N-sec/m ) 

N 

n 

number of orifice or capillary restrictors 

NH 


number of values of thrust clearance 

NW2 


number of values of ^2 

PS 

Ps 

2 2 

lubricant supply pressure, Ib/in. (N/m ) (needed only 
if CENT = . FALSE.) 

RC 


radius of restrictor circle, in. (m) 

RHO 

p 

lubricant density, Ib/ft^ (kg/m^) 

RI 


inner radius of bearing, in. (m) 

RO 


outer radius of bearing, in. (m) 

RPI 

V 

inner radius of pocket, in. (m) 

RPO 


outer radius of pocket, in. (m) 

TMAXX 


maximum bearing torque for which calculations are 
desired, in. -lb (N-m) 

TMIN 


minimum bearing torque for which calculations are 
desired, in. -lb (N-m) 

W1 

0^1 

speed of upper thrust face, rpm 

W22 


initial value of it> 2 , rpm 


Any number of cases may be run at one time. 


Program Output 

Figure 5 illustrates program output. All input data are printed at the beginning of 
each case. After the input data, the speed of the upper thrust face is printed, followed 
by the supply pressure if that was calculated by the program (CENT = . TRUE. ). The 
following line shows the speeds of the lower and upper thrust faces, 0^2 00 ^, and the 

calculated mean speed (eq. (7)). The next line gives the torque due to the journal 
bearings and the loads the journal bearings will carry at an eccentricity ratio of 0. 5. 
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The following line comprises column headings. These headings include (1) clearances h^, 
hj, and hp, (2) total bearing load W and torque T, (3) total lubricant flow Q, (4) load 
contributions of outer and inner thrust face lands, and Wj, (5) load contribution of the 
pocket for r > and r < R^, (6) torque contributions of the outer and inner thrust face 
lands and pocket, (7) lubricant flow rates outward and inward from orifice radius R^ (in- 
ward flow is negative), (8) pressure p^ downstream of the restrictors and pressures 
Pp. and p. at radii Rp^ and R^^, (9) pressures Pp^ and p^ at radii Rp^ and R^, 

(10) film rotational Reynolds numbers for outer and inner thrust face lands and pocket, 
and (11) radius R and pressure and the derivative dp/dri „ . If R^ < 

"min I ^ ^o ^min 

Rq or dp/dr > 0 the pressure somewhere in the bearing film is subambient. Under 

this condition, a real bearing will cavitate; computer results are inaccurate. Pressure 
Pmin found by extrapolating from r = R^ until dp/dr = 0 using equation (8). 


SAMPLE RESULTS 

The pressurized thrust bearing for which results are given is illustrated in figure 1 
as part of a series -hybrid rolling-element bearing. This bearing was evaluated experi- 
mentally in reference 1. The centrifugal field in the rotating upper thrust member was 
used to pressurize the lubricant. Bearing dimensions are given in U. S. Customary 
Units in the computer program output of figure 5. Figure 6 shows bearing thrust load as 
a function of outer thrust face clearance h^. Typical of externally pressurized bearings, 
the stiffness becomes very small as the clearance approaches zero. (Stiffness is given 
by the slope of the load -clearance curve. ) When the lower thrust face rotates at the same 
speed as the upper face, the load at zero clearance is somewhat higher than when the 
lower thrust face is stationary. At larger clearances, however, the load is less when 
both faces rotate. Centrifugal force is responsible for this difference. The rotative 
speed used to calculate centrifugal effects (^^) is between the upper and lower thrust 
face speeds (eq. (7)). Thus, with both faces rotating, centrifugal effects are greater 
than with only one rotating. 

Figure 7 gives the variation of lubricant throughflow with change in clearance. As 
expected, flow increases with increasing clearance. Increased centrifugal effects in- 
crease the slope of the curve when both thrust faces rotate. 

Figure 8 shows the torque of the bearing assembly as a function of clearance. 

Torque decreases with increasing clearance. When the lower thrust face rotates at half 
the upper thrust face speed, torque is halved compared to the case when 0 O 2 = 0, since 
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torque is directly proportional to the speed difference between the two faces. When the 
two faces rotate at the same speed, the torque is, of course, zero. 

Figure 9 shows the variation of load with lower thrust face speed for various bear- 
ing torques. This has been plotted from the alternate form of program output, in which 
clearance is varied to give preassigned torque values. This form is convenient for deter- 
mining bearing speed for given load and torque, as was necessary in reference 1. 

In figures 6 and 9 the curves stop before reaching zero load. This is because, as 
the bearing clearance increases, cavitation eventually sets in, making load calculations 
inaccurate. The presence of cavitation is indicated by negative intermediate pressures 
in the program output. 


Comparison with Experiment 

Figure 10 compares speed sharing predictions of the analysis with experimental data 
from reference 1. Oil viscosities corresponding to experimentally measured ball- 
bearing outer -race temperatures were used in the computer program. Points for the ana- 
lytical curve were obtained from plots similar to figure 9 using measured bearing torque. 

The experimental apparatus had a centrifugal oil supply. Therefore, supply pressure 
varied with shaft speed according to equation (18). When the supply pressure becomes 
high enough (high enough shaft speed) to lift off the fluid -film thrust bearing, the inter- 
mediate speed drops abruptly, according to the analysis. After lift-off, intermediate 
speed rises at a slightly lower rate than shaft speed. 

Agreement is generally good between analysis and experiment, though experimental 
lift-off speeds and intermediate speeds are somewhat higher than predicted. Possible 
reasons for this are (1) the analysis assumes a line source of lubricant in the fluid -film 
thrust bearing, whereas the experimental bearing has only four orifices, and (2) the fluid- 
film bearing may be cooler than the ball bearing outer race, resulting in higher oil vis- 
cosities than assumed. Fluid turbulence probably did not affect the results significantly, 
since the highest Reynolds number predicted by the analysis is 1020 for the data of figure 
10. This is only slightly above the threshold of 1000 reported in reference 6. 


CONCLUDING REMARKS 

An analysis has been presented and a computer program has been developed to en- 
able rapid evaluation of rotating pressurized thrust bearing designs using an incompres- 
sible lubricant. Included in the analysis are the effects of two self-acting journal bear- 
ings which may be used to provide a radial load capacity. Bearing load, torque, lubri- 
cant flow rate, and other quantities of interest are calculated. Either orifice or capillary 
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compensation can be used and effects of bearing rotation are included. Program input 
and output can be in either U.S. Customary or metric (SI) units. The computer program 
was written in FORTRAN IV; it can be used on most modern digital computers. 

Analytical predictions agree well with experimental, data from a series -hybrid fluid- 
film rolling-element bearing. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, October 12, 1970, 

129-03. 
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APPENDIX - COMPUTER PROGRAM FOR ANALYSIS OF 


PRESSURIZED ROTATING THRUST BEARINGS 

$ IBFTC RCFF 2 

rf:al li,lo,lr,mu 

LObICAL -)tLT,LENr,MFTRIC 

NAM EL I ST /BRb / CO ,C I tC 0 fJ ,DH ,DT ,OW2 tHI I , HOO, H PP , 1. 1 , LOt L«, MUt N, NH, 

1 Nw2,PS,RC,RHO,RHL;r'1,RI ,RO,RPl ,RPO, T MAXX , T Ml N, W I ,K2 2 . CENT , M ET R I C 
C CD = LiRIFICG DISCHARGE COEFFICIENT NH = NUMBER CF MG'S 

C Cl = INNER JNL bRC RAO CLtARANCE Nw2= NUMBER OF W2*S 

C CU = OUTER JNL ERG RAD CLEARANCE PS = SUPPLY PRESSURE, PS I (N/M21 

C D = ORIFICE OR CAPILLARY DIAMETER RC = RADIUS CF ORIFICE CIRCLE 

C DF = INCREMENT IN HU, HI, AND HP RH0= LUBRICANT DENSITY, LB/FT3 (KG/M3) 

C OT = increment in TORQUE RI = INSIDE BEARING RADIUS 

C DW2= increment IN W2 , RPM RO = OUTSIDE BFARING RADIUS 

C FI1= INITIAL INNER THRUST CLEARANCE RPI= INNER POCKET RADIUS 

C hOC)= INITIAL OUTER THRUST CLEARANCE RPO= OUTER POCKET RADIUS 

C HPP= INITIAL POCKET CLEARANCE TMAXX= MAXIMUM TORQUE, IN LB (N Ml 

C LI = inner JOURNAL BEARING LENGTH TMIN= MINIMUM TORQUE, IN LB (N M) 

C LU = LUTER Ji-JURNAL BEARING LENGTH = SHAFT SPEED, RPM 

C LR = CAPILLARY LENGTH W22= INITIAL LUwER THRUST FACE SPEED, RPM 

C MU = LUBRICANT DYNAMIC VISCOSITY CENT IF. TRUE. PS CALC FROM CENT FORCE 

C N = NUMBER JF ORIFICES METRIC IF. TRUE. USE METRIC UNITS 

c ALL Lengths ake inches (meters) 

F(E) = E/ A./ ( 1 .-c=!'E )<‘'*‘2*SQRT( PI*PI «( 1 . -E«e)-H6 . *£*F ) 

C F(t) IS FAuTUR Tu COMPUTE JOURNAL BRG LOAD (FROM SHURT 3RG THEORY) 

DATA P1/3.141bS26B/,CO/.6/,CO,CI/1. ,1./,HII , HOC, HP P/0 . , 0 . , . 1/ 

DATA Oh/. 00)1/, LI ,LQ,LR/0. ,0. ,0 . / , rt2 2 /O . / , MET R I C/ . F AL S E ./ 

OATA TMAXX,D T , TMI N/0. ,0. ,U. / 

I CM AX = 5 

C ICMAX = MAXIMUM NUMBER OF ITERATIONS TO GET TORQUE WANTED 

5 READ ( 5, BRG) 

DHP = FPP - HOO 
DHI = HI I- HOO 
EPST = .01 

IF (LR .E ) .0. ) (vRi T£ (6,6) 

IF (LR.NE.O.) WRITE (6,7) 

6 FORMAT (^SHl ANALYSIS OF ORIFICE COMPENSATED THRUST BEARING ) 

7 FORMAT (SlHl ANALYSIS OF CAPILLARY COMPENSATED THRUST BEARING ) 

R H J M = R H 0 

IF ( .NOT. METRIC ) RHOM = RH 0/1 72 B . /3 36 . 

RC2 = RC-XRC 
R0 2 = RC'-XRO 
RI2 = RI'f'RI 
RP I 2 = RP I=XRP I 
RP0 2 = RPCj^RPO 
WIR = Wl>:=PI/30. 

IF (CENT) PS = RHOM=i‘RC2*WlR*WlR/2. 

WR I TE ( 6, BRG ) 

IF (CENT. AND. .NOT. METRIC) WRITE (6,10) «1,PS,RC 
IF (CENT. AND. METRIC) WRITE (6,11) Wl.PS.RC 

IF (.NCT.CcNT) WRITE (6,12) Wl 

10 FORMAT (5HK Wl F 1 0 . 1 , 5H RPM, A-X 2 7HC E NTR I FUG AL SUPPLY PRESSURE 
1 G10.3,11HPSI AT RC = G10.3 ,6HI NCHES) 

11 FORMAT (SHK wl F10.1,5H RPM, 4X 27HCENTRI FUGAL SUPPLY PRESSURE 
1 G1C.3,13H iN/M2 AT RC = G10.3,7H METERS ) 

12 FORMAT ( 3HK Wl F'10.1,4H RPM ) 

RKO = CD*FLJA TIN )*PI*D*D/4.*SQRT(2. /RHOM) 
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SKU = HK.j«'<KU 

XC = 12f'.*MU=!'LH/FLUAT(N) /PI /n**4 
SIX = e.*MU/PI 
ROPL = ALOG(KO/f<PL)J 
KPJCL = AL0GIRP3/KC) 

KCP IL = aLJG(RC/RPI ) 

RP IL = AL'.IGIRP I /R I » 

C FACrUKS FCK CUMPUTI.MG THRUST 3EAP1RG LOADS AND LUURICANT FLOrfS 
PF;; = RHJM/^.*(RU2 - RP02 ) 

PF3 = KHUH/2.* (KP02 - RC2 ) 

PF4 = KHUH/2.* (RPI2 - RC2 ) 

PF5 = RHUM/2.*(RI2 - kP12) 

V-iOF = S1X*(K0 2*RUPL- ( RC2-RP02 » /2. ) 

WPJF = SI X<=(FPO2*kPO0L - (RGP2-RC2 )/2. » 
aPIF= Slx=i‘( (RC2-RPI2) /2. - RPI2*RCPIU 
WIF = SIX*((KPI2-KI2) /2. - RI2>!‘RPIL> 

W 2 = W 2 2 
F = . 5 

C FACTORS FOR UJ.'IPUTING JOURNAL BLARING TURQULS f LOADS , AND FLOW RESISTANCES 
Fb = F(LI 

TJuF = 2.*PI*ML*RC2<‘RC*LO/CO 
TJIF = 2. *PI *MU*RI 2*R 1 /C H‘L I 
mJIF = FL<=MU=t'R I /C I /Cl *L I 
WJOr = MU*kO>i‘LC*^‘3/CO/CO«FE 
C DO LOOP FOR KANt.L OF SPEEDS 
DO 50 Nvv= I,Ni^2 
HO = HOD 
HI = h II 
HP = FPP 

WU2 =rtl*W2 + . 3=!'( WI-W2) **2 

WO = SORT (WO 2) 

rtk ITE (6, 15) w2,Wl .WO 

15 FORNAI (5HL W2 F10.I.5H KPM, 3X 2HK1 F10.1.5H RPM, 3X 2HW0 
1 F 10. l.AH RP'1 ) 
rtOR = n02«PI>!'P I /9C0. 

W2k = k-2*PI/30. 

OWR = V>1K - W2R 
TJO = TJIF^^DWR 
TJ I = IJlF^JwR 
WJU = «JOF*l)wR 
W J I = *> J I F D A R 

IF ( .NCT.PbTR 1C ) WRITK (6,201 T JC , T J I , E , w JO , W J 1 
IF (NETRIC) WRITE (6,21) T JO, T J I , E , w JO, W J I 

20 FORMAT (a2H journal BEARING TORQUES TJO.TJI, IN LB 2G10.3, 

1 22H LOADS AT ECC RATIO = F5.2.14H WJO.WJI, LB 2G10.3 ) 

21 FORMAT (40H JOURNAL BEARING TORQUES TJO.TJI, N M 2G10.3, 

1 22H LOADS AT ECC RATIO = F5.2.19H WJO.WJI, NEWTONS 2G10.3 ) 

IF ( .N JT.METR IC ) WRITE (6,25) 

IF (MEIRIC ) WRITE (6,26) 

25 FORMAT ( 9HkH0,H1,HP 9H LOAD.LB 2X 12HFLUW, IN3/SEC 2X BHWO.WI.LB 

1 3X 10FWPU,WR1 ,Li3 3X 3HTU,TI,TP 5X 5HQ0,QI 5X 9HPC,PPI,PI 5X 

2 6HPP0,P.T 2X IIHKEYNILDS NR IX 12HR( PHI N) , PMI N / 

3 3X 4HMILS 4X 1 2H TUr. J UE , I N LB 36X 5HIN LB 6X 7HIN3/SEC 7X 3HPS I 

4 9X 3HPS1 4X 1 IHIEO ,RE I ,REP 3X 8H0P0R(R0) ) 

26 FORMAKllHK HU.Hl.HPdH LOAD.N 3X 12HFL0W, M3/SEC 3X 7HW0,WI,N 

1 4X VHWPJ,rtPl,N 4X 8HT0,TI,TP 5X 5HQ0,QI 5X 9HPC,PP1,PI 5X 

2 6HPP0,P0 2X IIHRFYNOLDS NR IX 1 2HR ( PMI N) , PMl N / 

3 6X 2HRM 8X I IH TORQUE, N M 36X 3HN M 7X 7H M3/ S EC 6X 5HN/M2 

4 7X 5HN/-12 3X 1 IHkc 0 , RE 1 , RE P 3X 8HD'»0R(R0) ) 
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PR3 = 

PF3^W3R 


PR2 = 

PF2<^WJR 


PR4 = 

PF4«W0R 


PR5 = 

PF5<=W0R 


PRO = 

PR2 + PK3 


PRI = 

PR4 + PR5 


TTUF 

= PI *MU=i'0WR/2. 


TTPF 

= TT0F^(RP02*KPG2 

- RPI2*RPI2) 

TTIF 

= TTGF«(RP12*RPI2 

- RI2^RI2) 

TTOF 

= TT0F^(R02=!'RG2 - 

RP02«RP02 > 

REF = 

OWR^RHOM/MU 



C FACTOR FOR COMPOTING FILM ROTATIONAL REYNOLDS NUMBERS 
DELT = .FALSE. 

TMAX = TMAXX 

C DO LOOP FOR RANGE OF CLEARANCES 
DO AO I=1,NH 
IC = C 

2 9 H0 3 = HO-^4=J 

HI3 = HI^AS 
HP3 = HP’)=<=3 
HIM3 = AMINl (HI3.HP3I 
HOM3 = AM INI (HO 3, HP 3) 

Xl= S I X*LO*HOM3/RO/CO*'i=3 
X2 = SIX^KOPL 

IF( HO 3 .NE .HOM3 ) X2 = X2<‘HCM3/H03 
X3 = SIX^=RPOCL 

IF( HP 3 .NE .HOM 3 ) X3 = X3*HOM3/HP3 


XA = SIX*RCPIL 

1F( HP 3 .NE .HIM3) XA = XA«^HIM3/HP3 
X5 = SIX^RPIL 

IF( HI 3 .NE .HIM3I X5 = X5*HIM3/HI3 
X6 = SIX *L1’«=HIM3/RI/CI**3 
XU = X I <- X2 + X3 
XI = X« + X5 + X6 
hXI = FIM3/XI 
HXO = HOM 3/XO 
HX = HXI + HXO 
IF (LR .EG .0. > GO TO 293 
CALCULATE PC FOR CAPILLARY COMPENSATION 

PC = (PS/XC - PRO=«^HXO - PRI*HXI ) / (HX l./XCI 
G = (PS- PCI /XC 
GO TO 2S5 

293 PHX = (PRU«=HXO + PRI*HXI)<'HX 
CALCULATc PC FOR ORIFICE COMPENSATION 

PC = (-PHX - SKO/2. + RKO»SQRT(PHX ^ SKU/A. «- HX*HX *PS I )/ HX/ HX 
Q = RKOASQRTIPS-PC) 

295 GO = ( PC + PRO ) /XO 
a I = - (PC PR 1 1 /XI 
PP I = PC + a IX=XA + PR A 
P IN= PPI + a I*X5 *■ PR5 
PPO= PC - QU*X3 + PR3 
PO = PPG - aO<‘X2 ♦ PR2 
WO = -ao*woF 

IF( HO 3 .NE .HUM3 ) WO = WO^^HUM I/HOO 

WO = PI*(WO F (PPO + PR2/2.) *(R02 - RP02I) 

WPO = -aO*WPOF 

IF( HP 3 .NE .HJM3 ) WPO= wPQ*H0M3/HP3 

WPO = PI*(WPO+ (PC + PR3/2.I *(RP02 - RC2 I I 

WP I = Ul^i^WPIF 

IF( HP 3 .NE .HIM3 ) WP I = WPI*HIM3/HP3 
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WPI= PI*(wPI + (PC + PK4/2.) RC2-RP[2n 
WI = Q 

IF ( HI 3 .NE .HIM 3) WI = Wl^HI M3/HI3 

WI = PI«(WI + (PPI 4- PK5/2.) ^(RPI2- RI2)) 

FAC = d.^MU^’^QO/PI 

IF( H02.NE .H0M3) FAC = FAC*H0M3/H03 

QO = Ua^HUM3 

0 I = Q I^H IM3 

TTO = 1.E6 

TTI = 1.F6 

TTP = 1.E6 

IF (HG.NE.O.) TTiJ = TTUF/HG 
IF (HI .NE.O.) TTI = TTIF/HI 
IF(HP.N£.C.> TTP = TTPF/HP 
T = TTO TTP + TTI TJO + TJI 
W = Wa + WP3 + WPI + WI 

CALCULATE MA)(IMUM FILM ROTATIONAL REYNOLDS NUMBERS FOR THRUST BEARING 
REU = REF^RJ^HO 
REI = REF^RPIX^HI 
REP = REF^RPO^HP 
OPDR = “FAC/RG + RHCM*^RG«WOR 
RM2 = FAC/RHOM/W'JR 
RM = SORT (AbS(RM2) ) 

PMIN = PPG - FAC=<=AHJa (RM/RPU) ^ RH OM^WOR/ 2 • ♦ ( R M2-R P02 ) 

IF (KM2.L T.O. ) RM = -RM 

IF (OELT. AND.AB S( T-TMAX) .GT.EPST.AND. W.GE.O. ) GO TO 35 
IF ( .NOT.METR IC ) WRITE (6,30) H 0 , W , Q , WO , W PG , T T 0 , QO , PC , P PO, R£0, RM t 
2 HI , T, WI, WP I , TTI ,QI t PPI , PG, RE I , PMI N, HP , TT P , P I N, R EP, DPDR 
IF (METRIC) WRITE (6,31) H 0 , W , Q , WO ,W PG ,T T 0, QO , PC , PPO, REO, RM, 

2 HI,T,Wl,wPI,TTI,QI,PPI,PG,REI,PMIN, HP , TT P , P I N , R EP , DPOR 

30 FORMAT ( /3PF8.3,0PG12.4,yG12.3/3PF11.2 ,0PG14,3, 9X 8G12.3/ 

2 3PFI1.2,47X 0PG12.3,12X G12.3,12X 2G12.3) 

31 FORMAT (/iPFlO. 5,0PG12.4,9G12.3/3PF13.4,0PG14.3,9X 8G12.3/ 

2 3PF13.4,A7X 0PG12.3,12X G12.3,12X 2G12.3) 

IF (W.LT.O.) GU TG 50 
35 HT = HO 

IF ( A3 S( r-TMA X) .LE .6PST.AN0. wl. NE. W2 ) GO TO 200 
IF ( T .GT. TMAXX • AND. . NOT.Dt LT. UR. WI. EQ. W2 ) GU TO 39 
IF ( IC .GT . ICMAX) GO TG ZOO 
iC = IC ^ 1 
DELT = .TRUE. 

IOC HO = HC ( TMAX-T) / (Tl-T)<^ (Hl-HO) 

IF (HO .LE .0. ) HO = . 005 
HI = HG + DHI 

HP = HQ + DHP 

IF ( ABS( r 1-TMAX) . LE. Ab 5( r-TMAX) ) GO TO 29 
HI = HT 

T1 = T 

GU TO 29 

20C TMAX = TMAX - D T 
IC = C 

DELT = .TRUE. 

IF ( TMAX. EQ. 0. .Ok . TMA X. LT. TMI N) GO TC 50 
GU TO 100 

39 HI = HU 
T1 = T 

HI = H I 4- OH 
HP = hP 4* DH 

40 HO = HG 4 - DM 
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bO W 2 = 2 

GO TO 5 
END 


4* DW2 
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ANALYSIS OF ORIFICE COMPENSATEO THRUST BEARING 


SBRG 


CD 

= 

6,0uC00C0E- J1 . 

: I 


l.OOOOOOOE-03, 

CO 

= 

1 .OOOOOOOE 

00, 

0 

= 

8. 9999999F-03, 

OH 

= 

1. 0000000 E-04 

OT 

= 

C. , 

JW2 


6.0000000E 03, 

HI I 

= 

4, 0000000 E- 

■02, 

HilO 

= 

0. , 

HPP 

= 

8.99999996-03 

L I 


S J1 , 

LO 


0. , 

LR 

= 

3 , 

» 

MU 

= 

9,4000noOE-07, 

N 


4 

NH 

= 

16, 


= 

3 , 

PS 

= 

1 .0227569E 

02, 

RC 

= 

1.2000000E 00, 

RHf! 

= 

6.0003003E 01 

R HUM 

= 

«.9R53942E-05, 

Kl 

- 

8. 7500000E-01 , 

R J 

= 

1 .6000000E 

00, 

RP I 

= 

9.249999ge-0I , 

RPO 

= 

1.37500006 00 

TMAXX 

= 

c. 

TMI N 


0. , 

ri 1 

= 

1 . 2000000E 

04, 

W22 

= 

0. 





CENT = 1, METK It = F , 


$ ENC 


Ul 12000.0 RPH, 

centrifugal supply 

pressure 102.3 PSI AT RC 

= 1.200 

INCHES 





V-2 0, RPM, 

Wl 12000, 

0 RPM, 

wO 6572. 

7 RPM 







JOURNAL BEARING TOKOUC S rJO,TJC 

, IN LB 

0 4 

.276 LOADS 

AT ECC RATIO 

= 0,50 

WJO,WJI, L8 0 

493. 

3 

1 2 

3 

4 

5 

6 

7 

8 

9 

10 


11 

HO, HI, HP LOAUtLO F LQ v, , I N 3 / SC C 

WO, Ml ,L0 

WPOtWPl , LB TO,T !,T P 

03, QI PC,PPI,PI 

PPO,PO REYNOLDS NB 

RI PMINl ,PHIN 

MILS T'DROUE,IN LB 



IN L8 

IN3/SEC 

PSI 

PSI REO 

,REI ,R6P 

DPDRtRO) 

C. 444.35 

0. 474E-01 

106. 6 

145.5 

J.IOOE 07 

0 

97.95 

107.6 

0 


4.343 

40. 00 O.lOOfc 

07 

23.69 

168 .4 

0.677E-02 

-0.4746-01 

35.47 

0. 1076 - 05 

0.^456 

04 

-455.3 

9,00 




0.506 


83.55 


0.149E 

04 

-434.2 

C.ICO 444.01 

0.473E-01 

1 06. 6 

145,4 

55.28 

0.447E-C3 

97.88 

107.5 

19,24 


4.342 

40.10 60.14 


23.87 

168.2 

C.675E-02 

-0.473E-01 

85,39 

0. 107F-05 

0. 4466 

04 

-454.8 

9.10 




0.580 


33.48 


0.150E 

04 

-433.9 

C.200 441.57 

0.506E-01 

1 06. 0 

144,6 

27.64 

0. 356E-C2 

97.34 

106.9 

38.48 


4.332 

40.20 32.49 


23.72 

167.3 

0 .673E-02 

-0.4706-01 

84.86 

0.1076-05 

0,4476 

04 

-451.3 

9.20 




0,573 


82.94 


0.1526 

04 

-431.7 

C.300 434,46 

0.580E-01 

104.5 

142.3 

18.43 

0, 119E-01 

95.79 

105.4 

57.72 


4.304 

40.30 23.28 


23.28 

164.4 

0. 6726-02 

-0.4616-01 

83.31 

0. 1076-05 

0.4406 

04 

-441.1 

9,30 




0.56 7 


81.39 


0.1546 

04 

-425,3 

C.4C0 419.00 

0. 716E-01 

lOl. 1 

137.4 

13.82 

0.273E-01 

92.40 

102,0 

76.96 


4.243 

40.40 18. b6 


22.33 

153.2 

0.6706-02 

-0.442E-01 

79.92 

0,1076-05 

0.4496 

04 

-419.1 

9.40 




0.56 1 


78.00 


0,1556 

04 

-411.3 

(.5C0 39C.79 

0.9136-01 

94.90 

128.5 

11.06 

0.506E-01 

86.20 

95.81 

96,20 


4.120 

40.50 15.09 


20.58 

146.6 

0.6686-02 

-0.4076-01 

73.73 

0.1076-05 

0.451E 

04 

-379.8 

9.5C 




0 .55 5 


71.81 


0.157E 

04 

-385.7 

C.6C0 347.52 

0.115 

85.40 

114.9 

9.213 

0.798E-01 

76,69 

86.30 

115.4 


3.946 

40.60 14.04 


17,39 

129.4 

0 .66 76-0 2 

-0.353E-01 

64.22 

0. 1076-05 

0.45ZE 

04 

-321.7 

9,60 




0.549 


62.31 


0.1596 

04 

-346.5 

C.7C0 292.23 

0. 140 

73.26 

97.48 

7.897 

0,111 

64.53 

74.14 

134.7 


3.700 

40.70 12.72 


14.45 

107.0 

0.665E-D2 

-O.204E-O1 

52.07 

P. 1076-05 

0.453E 

04 

-251. S 

9.70 




0 .544 


50.15 


0.1606 

04 

-296.4 




Figures. - 

Example of computer program output 
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0 .01 .02 .03 

Outer thrust face clearance, h^, mm 


Figure 7. - Variation of lubricant flow with clearance. 
Upper thrust face speed, 12 000 rpm. 



CO 

cn 


Thrust load, 


500- 

Lower thrust 



0 .5 1.0 1.5 

Outer thrust face clearance, hg, mils 


0 .01 .02 .03 
Outer thrust face clearance, h^, mm 

Figure 6. - Variation of load with clearance. Upper thrust 
face speed, 12 000 rpm. 



Lower thrust 




0 .01 .02 .03 

Outer thrust face clearance, Uq. mm 


Figure 8. - Torque of bearing assembly. Upper thrust face 
speed, 12 000 rpm. 



■g 


■g 

E 


24 000 1— 


20 000 U- 


16 000 ^ 


12000 ^ 


8 000 


4 000^ 


Thrust load, 
lb(N) 

O 100 (445) 

□ 200 (890) 

A 300 (1340) 

Solid data points indicate 

theoretical lift-off speeds 
Open data points indicate 
experimental results 

Intermediate 
speed equals 



4000 


8000 12 000 16 000 
Shaft speed, rpm 


20 000 


24 000 


Figure 10. - Comparison of theoretical and experimental speed sharing 
of the series hybrid bearing. 
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